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Abstract 

In a previous paper (Avotina et al. Phys. Rev. B 71, 115430 (2005)) we have shown that in 
principle it is possible to image the defect positions below a metal surface by means of a scanning 

i-C tunnelling microscope. The principle relies on the interference of electron waves scattered on the 

defects, which give rise to small but measurable conductance fluctuations. Whereas in that work 
the band structure was assumed to be free-electron like, here we investigate the effects of Fermi 
surface anisotropy. We demonstrate that the amplitude and period of the conductance oscillations 
are determined by the local geometry of the Fermi surface. The signal results from those points 
for which the electron velocity is directed along the vector connecting the point contact to the 

t^ \ defect. For a general Fermi surface geometry the position of the maximum amplitude of the 

conductance oscillations is not found for the tip directly above the defect. We have determined 
q ■ optimal conditions for determination of defect positions in metals with closed and open Fermi 

surfaces. 

PACS numbers: 73.23.-b,72.10.Fk 
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I. INTRODUCTION 



The interference of electron waves scattered by single defects results in an oscillatory 
dependence of the point contact conductance G (V) on the applied voltage V. This effect 
originates from quantum interference between the principal wave that is directly transmitted 
through the contact and the partial wave that is scattered by the contact and the defect or 
several defects. Such conductance oscillations have been observed in quantum point contacts 
[1-4] and investigated theoretically in the papers [3, 5-7]. 

In our previous paper [8] the oscillatory voltage dependence of the conductance of a tunnel 
point contact in the presence of a single point-like defect has been analyzed theoretically 
and it has been shown that this dependence can be used for the determination of defect 
positions below a metal surface by means of a scanning tunnelling microscope (STM). In 
the model of a spherical Fermi surface (FS) the amplitude of the conductance oscillations is 
maximal when the contact is placed directly above the defect. The oscillatory part of the 
conductance AG for this situation is proportional to 



AG(V) -cos 2z \ k 2 F + 



2meV\ 



where zq the depth of the defect and kp and m are the Fermi wave vector and effective mass 
of the electrons [8]. Materials with an almost spherical FS are most suitable for this model. 

In most metals the dispersion relation for the charge carriers is a complicated anisotropic 
function of the momentum. This leads to anisotropy of the various kinetic characteristics 
[9]. Particularly, as shown in Ref. [10], the current spreading may be strongly anisotropic in 
the vicinity of a point-contact. This effect influences the way the point-contact conductance 
depends on the position of the defect. For example, in the case of a Au(lll) surface the 
'necks' in the FS should cause a defect to be invisible when probed exactly from above. 

Qualitatively, the wave function of electrons injected by a point contact for arbitrary FS 
e (p) = Ef has been analyzed by A. Kosevich [10]. He noted that at large distances from 
the contact the electron wave function for a certain direction r is defined by those points on 
the FS for which the electron group velocity is parallel to r. Unless the entire FS is convex 
there are several such points. The amplitude of the wave function depends on the Gaussian 
curvature K in these points, which can be convex (K > 0) or concave (K < 0). The parts of 
the FS having different signs of curvature are separated by lines of K = (inflection lines). 
In general there is a continuous set of electron wave vectors for which K = 0. The electron 
flux in the directions having zero Gaussian curvature exceeds the flux in other directions 
[10]. 

Electron scattering by defects in metals with an arbitrary FS can be strongly anisotropic 
[9]. Generally, the wave function of the electrons scattered by the defect consists of several 
superimposed waves, which travel with different velocities. In the case of an open constant- 
energy surface there are directions along which the electrons can not move at all. Scattering 
events along those directions occur only if the electron is transferred to a different sheet of 
the FS [9]. 

In this paper we analyze the effect of anisotropy of the FS to the possibility of determi- 
nation of the position of a defect below a metal surface by use of a STM. We show that the 
amplitude and the period of the conductance oscillations are defined by the local geometry 
of the FS, namely by those points for which the electron group velocity is directed along 



the radius vector from the contact to the defect. General formulas for the wave function 
and point contact conductance are obtained in sections II and III. In Sec. IV the asymp- 
totic forms of the wave function and the point contact conductance for large distances of 
the defect from the contact are found. The general results are illustrated for two specific 
models of the FS: an ellipsoid (Sec. V) and a corrugated cylinder (Sec. VI). Using these 
models, for which analytical dependencies of the conductance on voltage and defect position 
can be found, we describe the manifestation of common features of FS geometries to the 
conductance oscillations: anisotropy of a convex part ('bellies'), changing of the curvature 
(inflection lines) and presence of open directions ('necks'). 

II. THE SCHRODINGER EQUATION FOR QUASIPARTICLES 



Let us consider as a model for our system a nontransparent interface located at z = 
separating two metal half-spaces, in which there is an orifice (contact) of radius R centered 
at the point r = 0. The potential barrier in the plane z = is taken to be a delta function, 

U(r) = Uf(p)8(z), (1) 

where p = (x,y) is a two dimensional vector in the plane of the interface, with r = (p,z). 
The function / ( p) — > oo in all points of the plane z = except in the contact, where 
/ (p) — 1. At the point r = r near the contact in the upper half-space, z > 0, a point- 
like defect is placed. The electron interaction with the defect is described by the potential 
D (r — ro) , which is confined to a small region with a characteristic radius td around the 
point r . 

It is known that one can obtain an effective Schrodinger equation for quasiparticles in 
a metal from the dispersion relation e (p) , (the band structure) by replacement of the 
quasimomentum p (below for short we write momentum) in the function e (p) with the 
momentum operator p = ^V [9]. Here we do not specify the specific form of the dependence 
e (p), except that it satisfies the general condition of point symmetry e (p) = e (— p). For 
simplicity we assume that FS has only one sheet; there is only one zone described by the 
function e (p) . In the reduced zone scheme a given momentum p identifies a single point 
within the first Brillouin zone. The wave function if) (r) satisfies the Schrodinger equation 
with an effective Hamiltonian e (p), 

e (p) if,(r) + [e-U (r) - eV (z)} $ (r) = D (r - r ) V (r) , (2) 

where U (r) is defined by Eq. (1), V (z) is the applied electrical potential, and e is the 
electron energy. 

We consider a large barrier potential U. In this case the amplitude t of the electron wave 
function passing through the barrier is 

*(*,*)« xu ' (3) 

where v 1 ^ and v v z ef are the z-components of the velocity v = de (p) /dp of incident electrons 
(in) and electrons specularly reflected by the barrier (ref), respectively. Under condition 
of specular reflection the energy e and the component of the momentum tangential to the 



interface, pt = (p x ,Py) at 2; = are conserved. The components of the electron momentum 
perpendicular to the interface, pf (pt,e) and p r f (pt,e) are related by the equations, 

e {pT, p?) = e {pf , P f) = e; p^ 11 = pf = Pt. (4) 

The velocities vf and v T z have the opposite sign 

v in • N < 0, v rcf • N > 0, (5) 

where N is a unit vector normal to the interface laying in the half-space of the electron wave 
under consideration (N = (0, 0, 1) for z > and N = (0, 0, —1) for z < 0). We will assume 
that the crystallographic axes in half-spaces z ^ are identical. In this case the momenta 
and velocities for electrons incident on the barrier and for those transmitted through the 
barrier are equal. 

In general Eq. (4) may have several solutions, i.e. several specularly reflected states may 
correspond to an incident state with momentum p™. Such reflection is called multichannel 
specular reflection [11]. Below we assume that there is only one reflected electron state. 

In the limit of a small probability of electron tunnelling through the barrier, \t\ ■€ 1 
the applied voltage drops entirely over the barrier and we take the electric potential to be a 
step function V (z) = V Q (—z) . The reference point of zero electron energy is the bottom 
of the conduction band in the upper half-space, z > 0. The conduction band in the lower 
half-space z < is shifted by a value eV. We also assume that the applied bias eV is much 
smaller than the Fermi energy and in solving the Schrodinger equation (2) we neglect the 
electric potential V (z) . Equation (2) can be solved by using perturbation theory with the 
small parameter \t\ <C 1 [12]. In the zeroth approximation in this parameter we have the 
problem of an impenetrable partition between two metal half-spaces. 

We start by solving for the wave function ^°> (e, p t ; r) for a tunnelling point contact of 

low transparency, \t\ <C 1 without defects (D (r — r ) = 0). The wave function i[) (r), in 
zeroth order in the parameter |t| ^C 1, satisfies the boundary condition tp^' (z — 0) = at 
the interface, 

^[, 0) (e, p 4 ; r) = e iptp/h (e ip ^ z/h - e ip ^ z/h ) . (6) 



Let us consider an electron wave exp (ipr/H) incident on the junction from the lower 
half-space z < 0, so that v z > 0. In this half-space to first approximation in the parameter 
t the solution ■?/>' ' (r) of the homogeneous Schrodinger equation can be written in the form 
[12]: 

^(°) ( r ) = ^ 0) (r) + ^ (r) ; z < 0, (7) 

where the second term, tp^ (r) ex t, describes the changes in the reflected wave as a result of 
transmission through the contact. The wave function transmitted into the half-space z > 
is proportional to the amplitude t, 

^(°) ( r ) = v9 (+) (r) ; z > 0. (8) 

The function tp^ ' (p, z) satisfies the condition of continuity and the condition of conservation 
of probability flow at z — 0. For small |t| these conditions reduce to 

¥><"> (p, 0) = y?« (p, 0) ; (9) 

£ e W* = /(p)v9 (+) (p)0 ). (10) 



In the absence of defects (D = 0) the solution of the Schrodinger equation is given by [8, 12], 



oo 



,(+; , „ - -, - - . / / d v't jpf n _ n '\ Xp'tP+pz +) ( e >Pt) z )/ h 



S + >(p.-~)=t J j^f F ^-Pt) e{t K t,}l " (") 



— oo 

where 

on 

e Kpt-p' t )p/ti 



F(p t -p' t )= j dp ; (12) 

-co 

pi (e, p' t ) are roots of the equation 

e{ti,p?>)=e(p), (13) 

corresponding to waves with velocities v z (e, p' t ) > and v z ~ (s, p' t ) < 0. 

Let D (r) be a spherically symmetric scattering potential for a point-like defect, with 
a range ro that is order of the Fermi wave length Xp (D -^ D (0), the maximal value of 
D, when |r — ro| ^> ro). For a point-like defect (td — ► 0), the right hand side in the 
Schrodinger equation (2) can be rewritten as -D (r — ro) ^ (ro) [16]. This makes it possible 
to find a solution to Eq.(2) by means of the Green function Gq (r', r; e) of the homogeneous 
equation (at D — 0). The wave function scattered from the defect, ip (r), can be expressed 
in terms of the wave function <^ + ) (r) transmitted into the upper metal half-space, 

* (r) = y,(+) (r) + ^ (r„) - J( /; ro) y (14) 

where 

J (r, r ) = y dv'D (r> - r ) G+ (r', r; e) . (15) 

Because the Green function has a singularity at r — > r' Eq. (14) is correct if the integral 
J (r, r ) (15) converges in the point r = r . 

To proceed with further calculations we assume that the scattering potential is small 
and use perturbation theory in the interaction with the defect. This implies that we take 
|J(r ,r )| <C 1. The wave function solution, that is linear in \t\ and with a contribution to 
first order in D, is 

i> (r) = v^ (+) (r) + v? (+) (r ) J dv'D (r' - r ) G+ (r, r'; e) . (16) 

The Green function Go (r, r'; e) in zeroth approximation in the parameters D and \t\ should 
be calculated from the wave functions (6), 

f dp ^ W W> (r')Y 

G+ (r, r'; e) = / — ^ ^- '—. (17) 

V ; J (2irh) 3 e-e(p)-i0 K } 

This is the Green function of the Schrodinger equation in the half-space z > with hard-wall 
boundary conditions at z — 0. Substituting the wave function ip (r), Eq. (6), for z > in 
Eq.(17) we find, 

j_ r dpt eiMp - P >)/r l+tp j +) z/H M zl/h _ ip (-) z , /h 



n+ ( r -'. F \ - _ ^ l p -*Pz 'z'/n _ -%p\ 'z'/n . > / 



For < z < z' one should make the replacements z <-> z' and p z ' <-> —p z in Eq. (18); p z 
are given by Eq. ( 13). 

The main contribution to the integral in Eq. (16) comes from a small region near the 
point r' = r . Far from the point r = r (|r — r | ^> Td) the solution (16) takes the form: 

V (r) = v? (+) (r) + gG+ (r, r ; e) ^ (r ) , (19) 

where 

g = Jdr'D(r'-r ) (20) 

is the constant of electron-impurity interaction. 
III. POINT-CONTACT CONDUCTANCE 



The electrical current I (V) can be evaluated from the electron wave functions, ip, of the 
system through [13], 

I (V) = -^-3 J dp! p Q (v z ) [n F (e) -n F (e + eV)\ . (21) 



Here 



ex 

I P = dxdy Re i^*v z ^) (22) 

— oo 

is the density of probability flow in the z direction for the momentum p, integrated over a 
plane z = xconst, n F (e) is the Fermi distribution function, and v is the velocity operator, 
v = de (p) /dp. For the definiteness we choose in Eq.(21) eV > 0. At low temperatures the 
tunnel current is due to those electrons in the half-space z < having an energy between 
the Fermi energy, e F , and e F + eV, because on the other side of the barrier z > only states 
with e ^ e-p are available. 

After performing the integration over a plane at z ^> z , where the wave function (19) 
can be used, we find the density flow (22) becomes, 



I p = \t(e,p t )\\ 3 hR*({v^) 2 ) i/(e)+ (23) 



\t(e,p t )\ 9* 2 R 4 Rc f dp' t U +) (p' t )e ip '^ h ( p - ip i+\ P ' t)zo /K_ p -i P i-H P ' t ) Zo /H 

h J (27rhf vi +) (pi) - v { z - ] ( P { 



dp" r i(p' t 'p Q +pi +) (p' t ')z () )/h_ 

(2nhy ' 



where the mean value at energy e is obtained by integration over the surfaces of constant 
energy in momentum space, 

r dSp 

(•••) £ = £P ^'l P . (24) 



r dSp 

J Ivl 



= p- 



scaled by the velocity, |v| 



J| . In Eq. (23), pi + ' (p t ) and pi (p t ) are given by Eq. (13) 

and v (e) is the electron density of states per unit volume. 

Taking into account Eq. (23) we can calculate the current- volt age characteristics I (V). 

The conductance G (V) is the first derivative of the current / (V) with V and in the limit 

T — > we obtain, 

81 

G{ y) = w = eV {£) (/p w • (25) 

After integration over S p (24), Eq.(25) should be expanded in the parameter eV/e-p <C 1. 



IV. ASYMPTOTICS OF THE WAVE FUNCTION AND THE CONDUCTANCE 



In this section we find the wave function at a large distances from the contact, r ^> Xp, 
and an asymptotic expression for the conductance in the limit of a large distance between 
the defect and the contact, tq ^> Af and a small contact radius, R <C Af, where Af is the 
characteristic electron Fermi wave length. For R — > the function F in Eq. (12) takes the 
form [8], 

F(p t -p' t )=nR 2 , (26) 

and the wave function (11) can be written as 



¥> 



(+) 



p, z) = t (e, p t ) %R 2 



d Pt _ e i{p' tP +vi + \p' t )z)/h 



{2nhY 



(27) 



Let us consider the integral 



A (r,e) 



dp't 
(2irhf 



a »r(pt,r) 



(28) 



where T (p t , r) is the phase accumulated over the path travelled by the electron between the 
contact and the point r, 

T(p t ,v) = ^(p t p + p^(p t )z). (29) 

This kind of integrals appear in the expressions for the wave function (Eqs. (27), and 
(18), (19)) and for the conductance (Eqs. (23), (25)). At a large distance, r ^> Af, the 
exponent under the integral in Eq. (28) is a rapidly oscillating function and the integral can 
be calculated by the stationary phase method (see, for example, Ref. [17]). The stationary 
phase points p t = pf are defined by the equation, 

or 



dp t 



0. 



Pi=Pi 



(at) 



(30) 



With Eq. (29) we find: 



p + z 



dp. 



(+) 



Pt) 



dp t 



,(+) 



Pt 



(st) 

Pt=p t 



v z (Pt) 



0. 



(31) 



pt=p 



(st) 



For r ^> Af the asymptotic value A as (r,e) of the integral (28) is given by 



A a 



cos-# 



r,e 



2Khr-s/\Ko 



cxp 




sgni^o 



Here 



r (5,r) = r(pi st) , : 



(32) 



(33) 



is the phase (29) in a point defined by Eq. (31), K(e,p) is the Gaussian curvature of the 
surface of constant energy e (p x ,P y ,Pz) = £, and cos-#(r) = z/r is the angle between the 
vector r and the z axis. At the stationary phase points the curvature K (e, p) can be 
written as 

v . _ 



K (e, n) 



(34) 



pt=p 



(st) 



where A 



ik 



-i- . is the algebraic adjunct of the element 

ik \") 



m 



ik 



(P) 



d 2 € 

dpidp k 



(35) 



of the inverse mass matrix m _1 [15]; rii are components of the unit vector n = r/r. Note 
that for an arbitrary FS TOjfc(p) in the point pt = P* depends on the direction of vector r. 
It is follows from Eq. (31) that the velocity at the stationary phase point p\ s is parallel to 
the radius- vector r = (p,z). 

If the curvature of the FS changes sign, Eq. (31) has more then one solution p t = p£ s s 
(s = 1, 2...). In that case the value of the integral (28) is replaced by a sum over all points 
Pts, which in the limit of large distances is, 



A i 



r,e) 



E 



A! 



r,e 



(36) 



with Af (r,e) given by Eq. (32) for each stationary phase point s. It may also occur that 
Eq. (31) does not have any solution for given directions of the vector r, and the electron 
cannot propagate along these directions. These two energy surface properties result in com- 
plicated patterns of the distribution of the modulus of the wave function: 1) For directions 
for which Eq. (31) has several solutions a quantum interference pattern of the electron waves 
with different velocities should be observed. 2) When Eq. (31) has no solution for the se- 
lected direction of the vector r classical motion in this direction is forbidden and the wave 
function is exponentially small. 

For large values r, r' ^> Af the asymptotic behaviors of the Green function, Eq. (18), and 
of the conductance, Eq. (25) at r ^> Af, can be found analogous to the evaluation of the 
integral in Eq. (28). The slowly varying functions of the momentum must be taken in the 
stationary phase point. For the partial wave scattered by the defect that is moving towards 
the interface p z (p t ) in Eq. (29) must be replaced by p z ~ (pA In this case the stationary 
points have a group velocity v^' directed from the point ro towards the contact. From the 
central symmetry of the FS, e (p) = e (— p), it follows that the two stationary phase points 
for the function T (jp t , r) are antiparallel, p(+)( st ) 



.p(-)(st) = p (st)_ 



Next, we derive an asymptotic expression for the wave function (19), with r, |r — ro| ^> Af, 



for a symmetric orientation of the FS with respect to the interface, so that v 



(+) 






(+) 
Pz = -p. 

the form 



(-) 



and m^ 



0, if % 7^ k. Under these conditions the wave function (19) takes 



i> (r) » y,« (r) + £-<p™ (r ) E 



1 



i,/.' ^ v ,« + , (p2) A r(p-p«.^-,) 



«i +, (p!5 



-Af (p-A,,|2-aol;e) , 



(37) 



with 



■Y 



(+) 



irf J2 A T 



r,e 



(38) 



In Eq. (37) the velocity i>i is taken in the stationary phase points p\ s — p t s correspond- 
ing to the directions of the vector with coordinates (p — p ,\z ± zq\). 

We have assumed for the Gaussian curvature, Eq. (34), that Kq ^ in the stationary 
phase points pf . For those points at which K$ = the integral (28) diverges. This means 
that the third derivative of the phase T (pt,r) (29) with respect to p t must be taken into 
account. In Sec. V this is done for a model FS having cylindrical symmetry with respect to 
an open direction. Here we only note that the amplitude of the wave function cp^ (r) (38) 
in a direction of zero Gaussian curvature is larger than for other directions, and decreases 
more slowly as compared to the ~ 1/r dependence of Eq.(32). This results in an enhanced 
current flow near the cone surface defined by the condition K = [10]. The same effect 
also appears in the second term of the wave function (37), that describes the scattered wave: 
the amplitude of this partial wave is maximal in the directions of zero Gaussian curvature. 

If the FS has a flat part, i.e. Eq.(31) holds at all points of a FS region of finite area Sq, 
the associated electron waves propagate in the metal without decrease of their amplitude 
[10]. For the flat part the dispersion relation e (p) can be presented as e (p) = vop, where 
vq = const is the electron velocity. For such FS the asymptotic value of the integral (28) is 



Afl s (r,e) 



vozSa ( icr 
2 exp 



|v | (2nhY 



ftv 



(39) 



When the distance between the contact and the defect is large, r ^> Af, we obtain the 
conductance of the tunnel junction, using Eq. (23) and the asymptotic expression for the 
wave function (37), 



Ct — Crfl 



/ 



V 



\ 



2nh 2 ( v 



4 +) ) } i/(e F ) sS 



J2 ^ Af (r ,5 F , eV) Im Af (r ,e F , eV) 



(40) 



/ 



where Gq is the zero-bias (eV — ► 0) conductance of the junction without defect , 



G = eVR*h(\tf) £F V>(e F )((vM) 



(41) 



£F 



In deriving Eq. (43) we have assumed that eV <C £ F and tq ^> A F . Therefore, all functions 
of the energy e in Eq. (43) can be taken at e = £ F , except for the phase To (e, r). When 
eV < e F , 

r (r ,e F + eV)^T (e F ) + ^eV, £-° ~ -^ (42) 

C£ F C£ F £p A F 

and when the product (eVy£ F )(r /A F ) ^> 1 clearly the conductance (40) is an oscillatory 
function of the voltage V. Note that if the inequality eV <C e F is not satisfied, the value of 
the conductance (41) as well as the amplitude of its oscillations depend on V. The periods of 
oscillations are defined by the energy dependence of the function T (e, r ) and they remain 
the same for any voltage. Below presenting formulas for the conductance in different cases we 
don't expand arguments of oscillatory functions in the parameter eV/e-p. The obtained results 
properly describe the total conductance at eV <C e F and also can be used for the analysis of 
periods of oscillations at eV < ep. In Eq.(40) the Af (r ,£F,eV) denotes the function (32), 
in which To = To (i"o, £ F + eV) . Equation (40) shows that for a tunnel junction of small size 
the amplitude and period of conductance oscillations depend on the local geometry FS in 
those points for which the velocity is directed along the vector r . 

In the case of a convex FS there is only one stationary phase point satisfying Eq. (31) 
which allows simplifying the expression for the conductance (40). The oscillating part of the 
conductance, AG, can be written as, 



AG 



G. 

v 



i +) ) 2 ) v (ep) (27r/i) 3 r\ 



vi +) (e, p*: 



Yl A ik n 0i n k 

i,k=x,y,z 



sin(2r (r ,e F + e\/)) 



(st) 

Pt=p 



where no = tcq/tq. The phase of the oscillations in the conductance, 2Fo, is determined by 
the phase that the electron accumulates along the trajectory from the contact to the defect 
and back 

2r (r ,5) = ^p( st )r , (44) 

where 2 |p( st )| is the chord connecting the two points on the surface of constant energy for 
which the velocities are antiparallel and aligned with the vector r$. 

If the direction from the contact to the defect coincides with a direction of electron 
velocities for a flat part of the FS we should use asymptotic expression (39) for the function 
Ag s (r ,e) when calculating the conductance (40). For this case the oscillating part of the 
conductance, AG, is given by, 

AG_ ngSy z0 ^ / 2 (e F + eV)r \ ^ 



Go /( (+)\ 2 \ I |2 / w »-n4 V ^l V 



((^ +) ) 2 ) |vo|V(£ F )(2^) 



Note that r /|v | = zq/v z q. In the case when the FS has a flat part, the amplitude of 
the conductance oscillations for those special directions does not depend on the distance 
between the contact and the defect. 

In the next sections we shall consider two models of anisotropic FSs. The first model, of an 
ellipsoidal FS, illustrates the main features of the conductance oscillations for metals with a 
convex FS, i.e. with positive Gaussian curvature. The FS of the second model has the shape 
of a corrugated cylinder, which illustrates the effects of sign inversion of the curvature and 

10 




y (X z units) 

FIG. 1: Gray-scale plot of the modulus of the wave function in the plane x = for an ellipsoidal 
FS. The shape of the FS is defined by the mass ratios m x /m z =l, m y /m z = 3, and the long axis 
of the ellipsoid is rotated by 7r/4 around the x— axis, away from the y— axis. The coordinates are 
measured in units of X z = hj \j2m zz e. The position of the defect is ro = (0, 5, 15). 

the presence of open directions of the FS. These two models allow us to obtain dependencies 
of the conductance on the applied voltage and on the defect position in analytical form. 

V. ELLIPSOIDAL FERMI SURFACE 



For an ellipsoidal FS the Schrodinger equation (2) can, in fact, be solved exactly in the 
limit R — > and the wave function (19) and the conductance (25) can be found for arbitrary 
distances between the contact and the defect. For this FS the dependence of the electron 
energy e on the momentum p is given by relation, 



e(p) 



1 E 

2 ^ 

i,k=x,y,z 



PkPi. 

m ik 



(46) 



were Pi are the components of the electron momentum p, 1/mjfc are constants representing 
the components of the inverse effective mass tensor m _1 . The tensor m -1 can be diagonalized 
to the form {m -1 } 



ik 



m~ l 5ik so that the momentum- space lengths ^phE^rrii correspond the 
semi-axes of the FS ellipsoid e (p) = ep. 

For the ellipsoidal FS in absence of a defect (zeroth approximation) the wave function 



■r 



(+) 



r) (27) can be obtained by integration of Eq. (27) over momentum, 



<P 



(+) 



where 



* (e> Pt 



To (e, r) 



zwR 2 (2s 



,3/2 e ir 



*r (r)) 



h 3 Vdet m 



r 3 




E 

i,k=x,y,z 



^ik n i n ki 



(47) 



(48) 
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and 



* (e, P< 



h 



iUra, 



\ 



2m zz e 



m, 



i,k=x,y 



PiPk 

m ik 



+ 



m, 




(49) 



t=x,y 



Here, \x ik are the elements of the mass tensor p, which is the inversive to the tensor m -1 
(pm l = I, with I the unitary tensor; fi ik = Aik / det m _1 ). The Green function (18) takes 
the form 



Gn 



2ttV26 



exp [t\/^l (r — r' 



8£m zz zz' 

h 2 



h 3 Vdet m~ 



T?Jr 



8em zz zz' 
h 2 



exp (zTp (r - r')) 
r (r-r') 



(50) 



Using Eqs. (47) and (50) we obtain the wave function (19) to first approximation in the 
strength of the impurity potential. The modulus of this wave function is illustrated in 
Fig. 1 for a plane normal to the interface passing through the impurity and the contact. 
The long axis of the ellipsoid is rotated by n/A around the x— axis, away from the y— axis. 
An interference pattern is visible of partial waves reflected from the impurity with those 
emanating from the contact that is influenced by the anisotropy of the electronic band 
structure. 

The conductance in the limit of low temperatures, T — > 0, is obtained from Eqs. (23), (25) 
by integration over all directions of the momentum p and integration over the space coor- 
dinate p in a plane z = const (z > zq), retaining only terms to first order in g (i.e. ignoring 
multiple scattering at the impurity site), 



G c[[ (V) = Gf < 1 



12ttVo ( 2£ f) 



3/2 



1 



1 - 



r o (£F)i*o) 



sin 2r (e F + eV, r 



th 5 y/rn^ det [m x ] T^ (e F , r 
1 



x 



(51) 



To (^f,i"o) 



cos2r (e F + eV,r ) 



with To (e, r) given by (44). 

The amplitude of the conductance oscillations is maximal when To is minimal. For a 
fixed depth zq this minimum occurs when the defect position p is in the point p 00 with 
respect to the point contact at r = 0, where 



_ TTlzz/ Wlzx 

Po ° " Z ° I m zz /m. 



■-!) 



(52) 



The minimal value of the phase then becomes, 



"oo = r (e F + eV, p 00 , zq) = -zq yj 2 (e F + eV) m z 



(53) 



The phase Too corresponds to the extremal value of the chord of the FS in the direction 
normal to the interface. Gq 11 in Eq. (51) is the conductance in the absence of a defect 
(9 = 0): 

(54) 






2e 2 J R 4 4 



97ih 3 U 2 v/detfm" 1 ]^; 
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FIG. 2: Dependence of the oscillatory part of the conductance, AG, as a function of the position 
of the defect p m the plane z = zq for the same shape and orientation of the ellipsoidal FS as in 
Fig. 1. The coordinates are measured in units A 2 f = h/\/2m zz EF and the defect sits at zq = 5. 

Fig. 2 shows a plot of the normalized conductance, G° n (0) /Gq 1 , Eq. (51), for the contact 
as a function of the position of the defect, p , in the limit of low voltage, V — ► 0. We find 
that G cY[ is an oscillatory function of the defect position that reflects the anisotropy of the 
FS and the oscillations are largest when the defect is placed in the position p 00 , denned by 
Eq. (52). 

For the ellipsoidal model FS the wave function and conductance have been obtained 
exactly, within the framework of the model. For large distances r and r they transform 
into the asymptotic expressions, Eqs. (37), (40). We do not present the asymptotic form 
explicitly but it agrees to within a term proportional to T^ 4 to the exact from, Eq.(51). In 
Fig. 3 we compare the results for the calculations of the conductance by using the exact 
(51) and asymptotic (40) expressions. The figure confirms that for relatively small distances 
(Fig. 3(a)) the asymptotic formula still qualitatively describes the conductance very well and 
that for larger distances (Fig. 3(b)) the two results are in a good agreement. The parameters 
for the FS in Figs. 2 and 3 are the same as those for Fig. 1. 
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FIG. 3: Comparison of the oscillating part of the conductance for an ellipsoidal FS calculated in 
point On (52) of the maximum amplitude bv using the asymptotic (dashed curve) and exact 

= 4 (a) and zq = 10 (b) in units of X z p. 



the point p (52) of the maximum amplitude by using the asymptotic (dashed curve) and exact 
(solid curve) formulas. The depths of the defect are zq 



VI. OPEN FERMI SURFACE 

The second model FS we want to discuss has the form of a corrugated cylinder (Fig. 4), 
which is open along the direction p u , 






71 IX 



e > Si, 



(55) 



where 2it/b is the size of the Brillouin zone, and m is an effective mass. We further impose 
that the momentum perpendicular to the symmetry axis of the FS remains finite, 



P_l_max _i Pi. _i P_Lmini 



where 



P±i 



y/2 



me: 



Pa 



iin (e) = y/2m(e-ei) 



(56) 



(57) 



are the maximal and minimal radii of the cylindrical surface, respectively. As a consequence 
of rotational symmetry the Gaussian curvature K (34) of the surface depends only on pj_. 
The central part of the surface ('belly') has a positive curvature K > while the ends near 
the Brillouin zone boundary ('necks') have negative curvature. In the direction perpendic- 
ular to the symmetry axis there are two partial waves propagating with different parallel 
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inflection lines 



FIG. 4: Model of an open FS. 

velocities, Vi and v 2 , belonging to the parts of FS having opposite sign of K. Rotating 
away from the perpendicular direction towards the axis the two solutions persist but the 
two corresponding points on the FS move closer together until they merge at the curve 
defined by K = 0, the inflection line. For directions beyond this angle (i.e. for 9 < 9 C in 
Fig. 4) no propagating wave solutions exist. On the inflection line a unique solution with 
velocity v c is found. For 9 > 9 C there are two stationary phase points p± s (s = 1,2) that 
satisfy Eq. (31), corresponding two different velocities v^ + ) (p± s ) directed along the radius 
vector r. The larger value, p^2, belongs to the belly of the FS (K > 0) and the smaller one, 
Pi.ii belongs to the neck (K < 0). At the inflection line of the surface we have, 



d 2 p u (e,p ± ) 



dp 2 . 



0, 



(58) 



P±=P±o 

which defines the value of perpendicular momentum p_|_ . From this condition we obtain 

RLO = VP-LmaxRLmin- (59) 

The cone inside of which no propagating states exist is defined by the condition 

b 



V± 



(P±m ax -P±mm)sgn(v„), 



= cot 1? c 

P±=P±o 

where the components of the velocity v u and v± are given by, 



(60) 



Vn 



e\b 



sin (p„6) ; 



v± 



Pi 



m 



(61) 



In spite of the simplicity of the model FS (55), the integrals in Eqs. (19), (25) cannot be 
evaluated analytically. We can only discuss the asymptotic behavior for r ^> Ap- Qualita- 
tively this result should also be valid for r > Af- For the directions that have two stationary 
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phase points, having opposite signs of the Gaussian curvature, Eq. (40) acquires the form 
/ 



G op = G op 



1 fifcos 2 ^(n ; y^ 



V 



(2irh) 3 h / (vi +) ) 2 \ v (e F ) r 2 -,-'=1,2 J 



:X 



l^ s) (e F , no) K^ (e F , n ) 



cos (rj, s) (e F + eF, r ) + | (1 - s)) sin (r<''> (e F + eV, r ) + | (1 - *'))) (62) 

where Gq P is the conductance of the contact without defect given by Eq. (41), I5 (e, r) = 

r(pt t) ,r),and^ ) (6,n) = K(pt t) ,e). 

The appearance of the conductance oscillations depends strongly on the orientation of 
the FS with respect to the interface. Below we will consider two specific orientations, having 
the axis of the FS either perpendicular or parallel to the interface. 



A. Direction of open FS perpendicular to the interface 

When the iso-energy surface is open along the contact axis z the components of the 
momenta in Eq. (55) are p± = a/p1"+~p| an d p u = Vz- in this case the conductance of the 
clean contact (without defect) becomes, 

_ ire 2 R 4 m 2 ejb 2 
G o ~ mu2 ■ (63) 

From Eq. (31) the stationary phase points for the iso-energy surface are, 
1 / 4 

P±s = 2 ( ^max + P±min ~ ^ C ° t2 $+ ( 64 ) 



The angle d = aiccos(z/r) is defined by the direction of the radius vector r. The Gaussian 
curvature K and the phase T in the points (64) are given by relations, 



4 S) (e, n) = h2{ ^f^ {vl^ + Pi m m " ^ cot 2 *) 2 - 4pi max pi min ; (65) 

T (^r) = ^ sV / - 2 + l/ 2 + ^arcsmy^— — . (66) 

The angle d in Eqs. (64), (65) is contained within the interval d c < 1? < tt/2, where the d c is 
given by Eq. (60). 

The modulus of the wave function is plotted in Fig. 5. For the calculation of the wave 
function, Eq. (37), we used formulas (65) for the curvature and (66) for the phase in the 
asymptotic expression for the integral A as (32). Although, strictly speaking Eq. (37) is 
not applicable in the vicinity of the contact and near the defect, nor inside the classically 
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FIG. 5: Gray-scale plot of the modulus of the wave function in the plane x = for a warped 
cylindrical FS having the open direction along the contact axis z. The coordinates are measured 
in units of An = h/y/2me. The parameters used in the model are £\/e = 0.9, b\/2ms = 5.2, and a 
defect sits at ro = (0, 13, 18) . 

inaccessible region, Fig. 5 illustrates the main features of this problem. One observes the 
interference of the two partial waves with different velocities, the existence of a forbidden 
cone, the anisotropy of the waves scattered by the defect, and the enhanced wave function 
amplitude near the edge of the forbidden cone. 

At the inflection lines, where K = and 9 is given by Eq. (60), the square root in Eq. (64) 
is equal to zero. For this direction two stationary phase points merge and the electron veloc- 
ity is directed along the cone of the classically forbidden region. The asymptotic expression 
for the conductance (62) diverges at these points, which implies that the third derivative of 
the phase, Eq. (29), with respect to p± must be taken into account. When the vector r 
connecting the point contact to the defect lies along the cone of the forbidden region the 
conductance oscillations have maximal amplitude and the conductance takes the form, 



fop 



u 



1 



Cg_ 

me i 



h A zl 



1/3 



V / P-LmaxP_Lmin (p±max ~ RLmin) 3 SHI I 2F 00 - 



5n 



where 



| _ |' , .. - , ._ 2Z [ VRLmaxRLmin / P_Lmax P±maxRLmi n 

1 oo — J- o lRi_o 5 z — -— h arcsm \ 

OH \RLmax-P±min V 2m£ l 



(67) 



e=EF+eV 



and C is a numerical constant, C ~ 1.97-10 4 . The energy dependencies of p± ma , x and p± y 
are given by Eq. (57). 
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FIG. 6: Dependence of the oscillatory part AG of the conductance, as a function of the lateral 
position of the defect p in the plane z = zq. The open direction of the FS is oriented perpendicular 
to the interface. The coordinates are measured in units of Aiip = h/yJ2me-p. The parameters used 
for the model FS (55) are ei/ep = 0.9, b^/2meF = 6, and a defect sits at a depth of zq = 11. 

Fig. 6 shows a plot of the oscillatory part AG of the conductance G op (0), Eq. (62), as a 
function of the lateral position of the defect p for a fixed distance zq from the interface. The 
oscillation pattern has a 'dead' region in the center, corresponding to defect positions inside 
the classically inaccessible part of the metal for electrons injected by the point contact. This 
region is defined by the cone (60) and its radius p 00 = z / cot ~& c depends on the depth of 
the defect under surface. The oscillations in the conductance are largest when the defect is 
placed at the edge of the cone p = p 00 . In this case the defect is positioned in a direction 
of velocity belonging to the inflection line of the FS and the electron flux in this direction 
is maximal. 



B. Direction of open FS parallel to the interface 



The second orientation we want to discuss is that with the FS (55) having its open 
direction (the axis) parallel to interface, with p_i_ = (p x ,p z ) and p u = p y . The existence 
of a classically inaccessible region for this geometry leads to a strongly anisotropic current 
density in the xy - plane. The expression for the point contact conductance without defect 



Gq P becomes, 



/nrop 



e 2 R 4 7r(2e-e 1 ] 
AhWU 2 



(68) 
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The expressions for the phase, Eq. (44), and the Gaussian curvature, Eq. (34), now read, 
if > (e, r) = i f v /(x^ + z^(2m(e-e 1 \ s )) + ^ arcsin ^/\f\ ? ( 69 ) 



and 




(s) e^si^e ( , 2cot 2 e V 8e cot 2 9 

Ky(e,n) = (-1) — W 1 + — — - - — m , (70) 



4 (e — £i\ s ) V \ inEib 2 J me\b 2 



respectively. For this geometry we use spherical coordinates, with 9 the angle between the 
vector r and the y-axis. The variables Ai 2 have been obtained from Eq. (31), 



ll 2cot 2 # . 1XS // 2cot 2 #\ 2 8ecot 2 #. 
^=2| l + ^F + ^ V( 1 + 7^-J ~^eW>- (7l) 

The first stationary phase point, Ai, corresponds to a positive Gaussian curvature and the 
second one, A 2 , to a negative curvature. For 9 = i? c , Eq.(60), they become equal, 



Al =A2= 2^T f^imax- \MmaxPiminJ . ( 72 ) 

and the curvature (70) vanishes. 

Figure 7, acquired by using Eqs. (37), (69) and (70), illustrates the interference of waves 
with different velocities the electrons emerging from the contact and the interference with 
the waves scattered by the defect. For this geometry the classically inaccessible region is 
found near the interface to both sides of the contact. 

Figure 8 shows a plot of the oscillatory part AG of the conductance G op (0) as a function 
of the lateral position of the defect p for a fixed distance z from the interface. In this case 
two 'dead' regions appear symmetrically with respect to the center of the oscillation pattern 
along the open direction of the FS. The center of the pattern corresponds to a defect sitting 
on the axis of the contact, p = 0, for which Ai = 0, A2 = 1, sin# = 1, and cosy? = 1. At 
this point Eq. (40) takes the form 



G op (V) = Gl p <M + 



op 1 -, , ^g 



2p±max(£F + eV)z 



2p±min(gF + elQV . 



+Zp_|_ m axP_L min COS 



h 

RLmax(£F + eV)z + RLmin^F + eV)z 



ft 

VII. DISCUSSION 



We have analyzed the oscillatory voltage dependence of the conductance of a tunnel 
junction in the presence of an elastic scattering center located inside the bulk for metals 
with an anisotropic FS. These oscillations result from electron waves being scattered by 
the defect and reflected back by the contact, interfering with electrons that are directly 
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FIG. 7: Gray-scale plot of the modulus of the wave function in the plane x = for the warped 
cylindrical Fermi surface with the open direction along the y-axis and parallel to the plane of 
interface. The coordinates are measured in units of Am. The Fermi surface parameters are £\/e = 
0.9, b\/2me = 5.2, and the defect position is ro = (0, 10, 25). 

transmitted through the contact. The introduction of anisotropic electron movement beams 
the following implication: several points on the FS may share the same direction of the 
group velocity vector v whereas other directions for v can be absent. Two non-spherical 
shapes for the FS have been investigated: the ellipsoid and the corrugated cylinder (open 
surface) . 

Contrary to the case of a spherical FS [8] in the ellipsoidal model (46) the center of the 
conductance oscillation pattern does not need to coincide with the actual position of the 
defect but is displaced over a vector p 00 = z (m zz /m zx ,m zz /m zy ). When the STM tip is 
placed at this point the oscillatory part of the conductance is given by, 



AG(V) oc cos (-zoy/2{e F 



eV) m z 



(74) 



The oscillation period depends on l/m zz , the component of the tensor of inversive mass 
(35) for motion in the z-direction, and on the depth zq of the defect. This allows us in 
principle to map out the positions of defects, as long as the shape of the FS is known. Apart 
from the period of the oscillations, there is also information in the amplitude. Since short 
periods should correspond to small amplitudes this may be used for a test of consistency. 
However, quantitatively the amplitude is also influenced by unknown factors such as the 
defect scattering efficiency. The ellipsoidal FS is exceptional in that the problem can be 
solved exactly. This allows us to compare the calculation with the asymptotic approximation, 
and this shows that the approximation works very well for distances larger than Xp. 
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FIG. 8: Dependence of the oscillatory part AG of the conductance, as a function of the lateral 
position of the defect p in the plane z = zq. The open direction of the FS is oriented parallel to 
the interface along the y— direction. The coordinates are measured in units of Amf- The parameters 
used for the model FS (55) are £i/£f = 0.9, b^/2rnep = 6, and a defect sits at a depth of zq = 11. 



In the case of the corrugated cylinder (55) the open necks cause cones with opening angle 
2i} c (denned by the inflection line of the FS) to be classically inaccessible. If the orientation 
is such that the open direction is orthogonal to the surface this will result in a 'dead' region 
with radius Zq/ cot d c (60) where no conductance fluctuations can be observed. Thus, by 
measuring the size of this dead region we directly obtain the position of the defect. The 
oscillation amplitude will be maximal at the border of the dead region, since the current 
density will be highest in the direction of the group velocity at the inflection line. In 
analogy to a hurricane the 'eye' is surrounded by a ring of intense currents. Such rings of 
high amplitude oscillations have already been reported very recently in experiments on Ag 
and Cu(lll) surfaces [18]. For our model FS, along this border the oscillating part of the 
conductance is, apart from a phase factor, described by, 



AG(VJ oc sin — — — 

Oh \ RLmax-P±i 



+ arcsin 



"imax 



Pi. maxP_l_ i 



2mei 



(75) 



where p± max and p± m j n are the maximal and minimal radii of the surface of constant energy in 
the direction perpendicular to the axis of the cylinder (57), E\ is the amplitude of corrugation 
of the FS, and 27r/6 is the size of the Brillouin zone. Again we find that the depth of the 
defect is determining the oscillation period, so that for given FS parameters this information 
can be exploited to investigate the structure of the metal below the surface. 
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If the open direction is parallel with the surface the highest amplitude will be found with 
the STM tip straight above the defect. For p = the conductance oscillations are described 
by Eq. (73). Clearly, the oscillation pattern is more complicated than that from Eq. 74, since 
there are contributions from the belly as well as from the neck parts of the FS, plus a sum 
frequency. For small necks the signal will be dominated by the oscillation due to the belly. 

Although the two models presented in this paper are still rather artificial, they provide 
insights that are quite valuable for experimental work in this field. The most prominent 
conclusion is that the regular oscillations due to convex parts of the FS, that behave as for 
the isotropic FS discussed previously [8], will often be dominated by signals due to special 
directions. On any surface that features regions of zero curvature, the strongest conductance 
fluctuations will come from electrons travelling with the group velocity of that region. Not 
only does this hold for the inflection lines around the necks in the (111) direction of e.g. Cu, 
Ag or Au, it is also true for the almost flat facets in the (110) direction of the same metals. 
In the case of an inflection line the signal will decay as r -5 / 3 , whereas for a flat facet the 
signal does not decay at all. This effect can be exploited for imaging defects up to much 
larger depths than previously estimated. The particular shape of the FS for nearly all metals 
contain many detailed features that will allow us to check the validity of the conclusions 
drawn from the measured data. 
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